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В статті розглянуті підходи для підвищення інформативності оцінки ритму серця людини, 
введені деякі допущення щодо структури випадкового процесу, який визначає можливості 1 межі 
застосування математичного апарату до аналізу кардіоритму. Пропонується використовувати підхід 
до аналізу серцевої діяльності, заснований на розбитті масивів ритмограм на ділянки, відповідні 
певному функціональному стану людини. Подається алгоритм визначення детермінізму масивів 
ритмограм серця людини у вигляді детермінованого чи ймовірнісного процесів. 
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Вступ 

Аналіз серцевого ритму (СР) знаходить широке застосування при рішенні задач 
діагностики і прогнозуванні в клініці при визначенні  психофізіологічних 
характеристик людини-оператора в ергономіці, під час тренувального процесу в 
спорті і для швидкої психофізіологічної оцінки функціональних станів людини для 
надійного і оперативного контролю |1,2). Різні методи дослідження СР можна 
класифікувати таким чином. Перша група включає алгоритми, які розглядають СР, 
як випадкову величину. В цій групі вивчаються статичні параметри СР: середнє 
значення, середнє квадратичне відхилення, асиметрія, ексцес, а також різні похідні 
цих параметрів. Крім цього у психофізіології використовуються деякі спеціальні 
показники: індекс напруги, критерій 5АВЗ5, показники напруженості або організація 
гістограми (2-4). Для більшості цих методик при обробці СР характерна мала 
трудомісткість, що дозволяє проводити обчислення на мікропроцесорах, достатньо 
висока наочність, забезпечуючи легке розуміння і трактування отриманих 
результатів. Їх основним недоліком є низька інформативність, викликана відмовою 
від вивчення тимчасової динаміки ритмограми. Крім цього, як буде показаний 
нижче, можливі випадки, коли дві або декілька ритмограм, різні структури, що 
мають принципово, можуть відповідати одній і тій же функції розподілу, а, отже, не 
розпізнаватися на рівні випадкових величин. Більш ефективними є методики, які ми 
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відносимо до другої групи, які аналізують СР, як випадковий процес. При такому 
підході виникає необхідність у введенні деяких допущень щодо структури 
випадкового процесу, визначаючих можливості і межі застосування математичного 
апарату. При цьому, чим більш суворі ці обмеження, тим більше дієві алгоритми 
аналізу можна використовувати. 

Постановка проблеми 

Ми пропонуємо використовувати підхід аналізу серцевої діяльності, 
заснований на розбитті масивів ритмограм на ділянки, відповідні певному 
функціональному стану людини. Після цього ці ділянки ритмограм подаються як 
суми детермінованого і стаціонарного  ергодичних процесів. Виходячи з 
літературних даних (|5| і аналізу експериментальних результатів, можна зробити 
допущення про наявність гармонійних складових структури СР. Тому, в якості 
детермінованого процесу ми пропонуємо застосувати суми вигляду: 

т 


У(:)- ав ї У а; є 5іп( а. - (ф.)) 
. і-1 . 

де ас - математичне очікування детермінованого процесу, «аг, б; та Ф; - 
амплітуда, частота і фаза і-той гармоніки відповідно, п - число існуючих гармонік. 
Стаціонарний процес передбачається з нульовим математичним очікуванням 1 
обмеженою дисперсією. 

Кінцевий етап пропонованого нами методу аналізу ритмограмм полягає в 
знаходженні математичного очікування, амплітуд і частот гармоніки, для 
детермінованого процесу (фази містять украй мало інформації і мають тільки 
допоміжне значення при обчислювальних операціях). Слід зазначити, що у 
психофізіологічних дослідженнях достовірно вдається звичайно описати лише одну 
гармоніку, що має найбільшу амплітуду |6,7|. Інші ж фактично пригнічуються 
шумом. Таким чином із запропонованої нами структури опису СР, видно, що в 
результаті аналізу експериментатор одержує всього чотири різні показники: 
математичне очікування, амплітуду і частоту | основної гармоніки для 
детермінованого процесу, а для стаціонарного - дисперсію. Кінцевий етап 
пропонованого нами методу аналізу ритмограмм полягає в знаходженні 
математичного очікування, амплітуд і частот гармонік для детермінованого процесу, 
а також дисперсію для стаціонарного. Фази містять украй мало інформації і мають 
тільки допоміжне значення при обчислювальних операціях. 

Запропоноване розкладання СР на суму детермінованого і стаціонарного 
Гауссових процесів дозволяє також одержувати чіткі оцінки структури ритмограм, 
сприяє більш строгому визначенню залежності структури послідовності В-К 
інтервалів ЕКГ від всіляких чинників, що впливають на досліджуваного. Вже при 
вивченні властивостей серцевого ритму на достатньо абстрактному рівні вдалося 
знайти обгрунтування окремих особливостей його структури. 

Такий підхід до розкладання СР, забезпечує просту розробку програми-імітатора 
ритмограм, істотно полегшуючу відбір алгоритмів аналізу СР відшуканням областей їх 
застосування, а також дозволяє визначити ряд інших властивостей, що представляють 
інтерес для дослідника. Так, вивчаючи коефіцієнт ексцесу масивів СР, ми знайшли що 
цей показник стабілізується при довжині ритмограми біля 400 В-К інтервалів. 

Основна частина 

При побудові гістограм К-К інтервалів ЕКГ для визначення функціональних 
станів  людини-оператора ми часто знаходили наявність бімодальних або 
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полімодальних структур, природу яких, використовуючи відомі методи аналізу 
серцевого ритму, задовільно пояснити не є можливим |5,7|. При використанні вище 
запропонованого методу аналізу серцевого ритму з'ясовано, що функції розподілу 
СР визначаються відношенням квадрату амплітуди гармонійної складової до 
дисперсії білого шуму. При цьому, чим це відношення більше, тим більше виражена 
бімодальність, при відношенні близькому до 2 виявляється ефект, існуючий на 
практиці полімодальності. При зменшенні цього відношення чітко виявляється 
одномодальна структура. 

Аналізуючи ритм серця людини-оператора в біотехнічних системах, виникає 
необхідність в оперативному і коректному описі основних властивостей 
кардіоінтервалограми. При цьому слід враховувати, що психофізіологічний стан 
людини часто досить змінний, має значну варіабельність, а разом з ним змінюється 1 
структура ритму серцевих скорочень (6,3). 

Враховуючи вище викладене, виникла необхідність в розробці математичного 
підходу, що дозволяє з високою точністю проводити аналіз досить коротких, 
порядка 100 значень, кардіоінтервалограм. Виходячи з цього, був розроблений 
алгоритм визначення параметрів серцевого ритму, що забезпечує отримання 
асимптотично ефективних по відношенню до сімейства | функцій втрат, 
асимптотично нормальних і сильно спроможних оцінок. Крім цього, при достатньо 
загальних допущеннях, отримані оцінки сходяться з вірогідністю «1» до істинних 
значень. Більш того, оцінка частоти спектральної складової сходиться швидше, ніж 
функція 1/М, де М - довжина послідовності, що аналізуються. 

Запропонований алгоритм розроблений на підставі методу максимальної 
правдоподібності з при припущенні, що  кардіоінтервалограма | може бути 
представлена, як випадковий процес, математичне очікування якого змінюється в 
часі по гармонійному закону, а шум має нормальний розподіл (для збіжності по 
ймовірності достатньо незалежності шуму і кінцівки його дисперсії). Алгоритм 


полягає в наступному: на підставі отриманих функцій від 4 (числом б), де А - 


частота спектральної складової, будують функціонал Е(Л), Мінімальне значення 
функціонала по відповідає оцінці частоти. 


-вез 1с05"і- 1-205 1 С05) 4)" 
- Уа зіп) А- со5)- 2-355 "шзіп) А. созі "А 
и - вані і-А4- па зіп) - 4). 
Аз - Ху хісзіп)- ач " шхі)з-азіві" 2 
я з-- Я кі-єнаі-2 3 й хід -:С05/- 4 


в й хі? М 2 
Азз - -303. 1х)) 

їі Йо Аз 
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Е(2) - Яіз Я 33 
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Де (ху у - послідовність інтервалів часу між скороченнями серця. На основі 


знайденої оцінки частоти гармонійної складової, вирішуючи нижченаведену систему 
лінійних рівнянь відносно азіп(ф)і Бсо5(ф) 


зоре ЗА засо5(ф) - Аз - 0 
12 45іп(Ф) -- А; асо5(ф) -- А;з - 0 


Отримаємо оцінки амплітуди і фази гармонійної складової. На підставі 

отриманих оцінок, по формулі 
п о хУМа(х; - авіш()2-- ф)) 

знаходимо оцінку математичного очікування - |. Обчисливши квадрат різниці 
між отриманим гармонійним процесом 1 вихідною послідовністю, отримаємо оцінку 
шумової складової. 

Пропонований алгоритм реалізований на мові СІН-. Для вивчення збіжності 
пропонованих оцінок до істинних значень параметрів випадкового процесу залежно 
від довжини вибірки, а також для обліку можливої обчислювальної помилки було 
проведене більше 500 обчислювальних експериментів. Отримані результати показали: 

При роботі на комп'ютері достатньо користуватися змінними СР стандартної 
довжини. Час аналізу вибірки, яка складається з 100 значень не перевищує 12 сек. 

Помилка оцінки частоти дихальної складової кардіоінтервалограми (КІГ) не 
перевищує 590 з вірогідністю 0.95 при довжині вибірки 25 і більш значень. Помилка 
оцінки амплітуди дихальної складової не перевищує 1590 за тих же умов. 

Дисперсія пропонованих оцінок 1, а, ф, в середньому на 3090 менше ніж при 
використанні оцінок, отриманих методами спектрального аналізу. 

При визначенності основних підходів до оцінювання варіабельності ритму 
серця людини, ми зазначили, що методи оцінки часової області опираються на 
статистичні методи і направлені на дослідження загальної варіабельності серцевого 
ритму (ВСР). Методи оцінки частотної області - досліджують періодичні складові 
ВСР та інтегральні показники ВСР. До цих методів відносять автокореляційний 
аналіз та кореляційну ритмограму. Статистичні методи полягають у вимірюванні 
інтервалів, а також у порівнянні показників електрокардіограми. Вони дають 
кількісну оцінку варіабельності кардіоінтервалограми (КІГ), котра являє собою 
сукупність КВК-інтервалів ЕКГ, котрі відображаються один за одним. Жоден з 
підходів не є універсальним |б6,7. 

Вибір способу оцінки КП залежить від мети, яку поставив перед собою 
дослідник. Це може бути дослідження коротких та довгих періодів реєстрації ЕКГ, 
дослідження ЕКГ здорової людини або із патологією, дослідження ритму у різних 
функціональних станах людини або у стані фізичного і емоційного стресу. Не кожен 
з вище зазначених методів здатен працювати у всіх режимах. Ми пропонуємо підхід, 
який враховує всі попередні методи і який дозволяє якісно оцінити організацію 
роботи ритму серця (А) за правилом Г. Фкуеорстера для розрахунку ступеня 
визначеності (організації) функціонування систем проходить |9,10) у вигляді: 
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п 
Ва1-Нх /Нах, де Нтах - Іоб п, Н(х)-- й ює рі, 
ї«і1 


де Нтах - максимальна невизначеність системи, НИХ) - поточна невизначеність 
системи; п - кількість станів системи, Рі- ймовірність прийняття системою і - того 
стану. При такому підході якість роботи серця визначалася за допомогою наступної 
математичної формули |10). Величина А при цьому виглядає для визначеного 
масиву кардіоциклів наступним чином: 


пт; 
лай Хеа100!! 
ін16 


Суть цього методу полягає у визначенні кількості класів (16) у заданому масиві 
КВ-інтервалів (100-200), знаходженні кількості ударів серця у кожному з класів (пт і) 
та подальшому розрахунку змінної А, яка буде коливатися в діапазоні від 0 до І. 
Чим більше графік розподілу випадкової величини ВВ-інтервалів у класах нагадує 
рівномірний розподіл 1 має ймовірнісний вигляд, тим ближче значення А до нуля 1 
тим гірша організація серця на цьому проміжку. І, навпаки, чим більше розподіл ВВ- 
інтервалів нагадує нормальний розподіл 1 чим організація ритму серця ближче до І, 
тим роботи серця більше детермінована. Якщо всі ВВ-інтервали знаходяться в 
одному класі, то А- І. 

У програмі реалізовано два підходи до дослідження змінної А. В обох 
варіантах необхідно задати інтервальний час, за який бажано вивести значення А. 
Наприклад, ми задаємо період визначення масиву даних ВВ-інтервалів кожні 200 
секунд. Перший спосіб полягає у накопиченні даних, тобто у змінній А будуть 
зберігатися дані за весь проміжок часу дослідження. При другому способі у змінній 
А зберігатимуться лише дані досліджуваної області, без попередніх. Обидва способи 
реалізуються як графічно, так і виводом у консоль поточних значень цих змінних. На 
графіку виводяться різними кольорами ці два способи, а в консолі змінні виводяться 
у вигляді Ат з Хі; А2 з Х.; При стабільній роботі серця буде спостерігатися 
циклічність даних, що виводяться, тобто, наприклад, кожні 3 хвилини параметри А! 
та А2 будуть однакові, але кожні 200 секунд - різні. Тому дана програма ще корисна 
тим, що може допомогти дізнатися цей показник. 

Для реалізації програмного продукту була обрана мова програмування С-нЕ, 
вона є високоефективною та швидкою у виконанні заданих операцій. Для 
створювання графічних проектів було підключено модуль ОрепСї,. Для звичайної 
роботи з «ехе-застосунком можна не встановлювати цей модуль, він потрібен для 
редагування і написання проекту. ОрепСі, - це графічний пакет, а не математичний. 
Це означає, що усі побудовані графіки та моделі робились власноручно та 
попіксельно. В ОрепСї, немає вбудованих математичних функцій побудови графіків, 
як це є в Маріє чи Майаб. Програма перебуває у графічному циклі. Це зроблено для 
того, щоб можна було повертатися, переглядати та аналізувати результати. Елементи 
програми оцінки ВСР людини подані нижче. 
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М з 213 А - 9.324392 А? - 9.324392 м - 675.117 
М «427 А - 9.337391 А? - 9.09226939 м - 675.761 
М «638 А - 9.427997 А? - 9.398472 м - 680.329 
М з 832 А - 9.265255  А2? - 4.273679 м - 694.231 
М а 1929 А - 9.297141 А? - Ф9.216877 м - 792.332 
М з 1231 А - 9.178545 А? - Ф9.194368 мМ - 795.362 
М а 1427 А - 9.248155  А2 - 9.249352 м - 710.687 
М - 1629 А - 9.261316  А2 - 9.267894 м - 711.817 
М е« 1831 А - 9.299459  А2? - 9.165515 М - 712.561 


5схасі5зсісаї! тесНіноайз: 
НВ - 81.5879906 ЬБс5/тіп 
ТВА - 738.4091672 т5 

О - 2219.5383309 т572 
50мМммМ - 47.111977 т5 
8ВМ550 - 19.513371 т5 
СУ - 6.380264 з 

мм590 - 48 Бс5 

рмм50 - 2.485759 з 


Сеотесгісаї| тегсНпоаз: 

Мо - 700.000000 ьс5 

дАто - 36.697723 зб 

МВ - 332 т5 

мхОомп - 1.539837 

5І - 78.953796 я5 

ІМВ - 1190.535309 5/5 
МРІЕ - 0.004303 1/(ьс5ж5) 
РАРВ - 0.052425 з/Ьс5 


п КкТіаз55ом - 7 
М - 1931 

А - 0.279372 

М - 714.311584 
Моаагг - 7089 


Етапи роботи виведені на блок-схему роботи ПЗ (рис.1). На початку необхідно 
запустити ПЗ на виконання. Далі програма з вхідного файлу «ат.їхі» зчитує ряд 
кардіоінтервалів, визначеної у певної людини, щоб зберегти результати в базі. 
Досліднику пропонується обрати один з декількох режимів роботи програми, які 
будують графічні та геометричні методи. Програма буде працювати до того 
моменту, поки користувач не натисне кнопки дій, після виконання яких програма 
завершує свою роботу. 

Програма на заданих даних з вхідного файлу будує графік, який показує зміну 
тривалості кожної кардіоінтервалограми в мілісекундах. Програма надає можливість 
за бажанням додавати або зменшувати яскравість розмірної сітки для покращення 
аналізу, натискаючи «Р» або «-» на клавіатурі. Також для зручності було поділено 
запис інтервалів по 250 секунд. Графік можна прокручувати вперед. 

У режимі «побудови графіку зміни ЧСС» програма будує для заданого масиву 
даних з файлу графік зміни ЧСС на всьому проміжку запису. У даному режимі також 
доступна можливість збільшувати або зменшувати масштаб графіка для його 
кращого вивчення. 

У режимі «побудови гістограми» програма будує гістограму аналізу ВСР. Цей 
метод корисний в дослідженнях, де класи гістограми розбиті із кроком 50ме, за 
вимогами та правилами оцінки організації. 
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Запуск 
програми 


зчитування ЕЕ 
інтервалів 
з файлу та введення 
йма"зи кизанцієнста 


обрання режиму роботи 
ІЗ 


ІТобудова 
інтервал окріємім 


їМобудова графіку 
змийнихя ЧІССЄСС 


Побудова 
кістоюраічмом 


Побудова 
скатерограми 


ІТобудова графіку 
організації серякя 


Вивід у консоль 


| основних змінних 


Створення ЕНКГІГМІ, 
та заповнення 
бази 


Еінецщь роботи 
програми 


Рис.1. Блок-схема роботи програмного забезпечення для визначення організації 
ритму серця людини. 


У режимі «побудова скатерограми» програма будує кореляційну ритмограму 
(скатерограму), один з найточніших та інформативних геометричних методів аналізу ВСР. 

Проводиться вивід у консоль основних змінних ЕКГ, що необхідно для 
постановки діагнозу. Для цього необхідно подивитися на отримані статичні змінні, 
які виводяться у консолі. Програма виводить такі змінні, як: НЕ, тавВВ, 0, 50ММ, 
ВМ55Ю, СУ, ММ50, рмММ50, Мо, Ато, УВ, МхОмМлп, 51, ТМВ, УРЕ, РАРК та інші. Всі 
необхідні дані виводяться у вигляді НТМІ,таблиці. 

На даному етапі 4 програма перебуває | у зд режимі тестування та 
експерементального дослідження для знаходження нових закономірностей. 

Висновки 

Враховуючи наведене вище, можна зробити такі висновки: 

1. Розроблено алгоритм визначення параметрів серцевого ритму, що забезпечує 
отримання  асимптотично ефективних, асимптотично нормальних |і сильно 
спроможних оцінок масивів ритму серця людини. 
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2. На підставі методу максимальної правдоподібності, розроблено алгоритм 
при припущенні, що кардіоінтервалограма може бути представлена, як випадковий 
процес, математичне очікування якого змінюється в часі по гармонійному закону, а 
шум має нормальний розподіл. 

3. Попередні дослідження показали, що при поточній оцінці масивів серцевого 
ритму достатньо користуватися змінними стандартної довжини. Час аналізу вибірки, 
яка має до 200 значень не перевищує 12 сек. 

4. Розроблено алгоритм та програму визначення ступеня детермінізму 
кардіоритму при різних психофізіологічних умовах, що дозволяє отримати значення 
організації отриманих масивів ритму серця людини. 
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ФуЛлСве Мао 


КЕ5ОМЕ 

Уи.Р. Согео, Н.А. ОтеїапепКо 

Меуу  аїсогійпт5 ої апа!узія тап'8 сагдіас гпуїпт 

Іп Фе агісіе Шеге аге Ше соп5ідегед арргоасрез Гог Бе гі5е ої іпіогтіпє ОЇ 
е5ійтацоп ої груїбта ої плап'з Беагі, їпе епіегед 5оте аззитріїопя їп геіайоп їо Фе 
зігисіите ої ассідетка! ргосе55 міс Фегегтіпез ро5516ійПез апа зсорез ої арріїсайоп ої 
таїфетайса! уерісіе їо Ше апаіузія ої сагфіогуШтя. І 15 5зибєезіед іо (аке арргоасі іо 
Фе апаЇузіз ої сагдїас асіуїсу, Баз5ед оп Іауїпе ої аггауб ої сагйфіоєгат оці оп агеа5 ргорег 
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го Фе дейпіге Гипспопа! 5(аге ої плап. ТБе аїсогійт ої дегегтіпацоп ої деїегтіпізт ОЇ 
аггауз ої груФт єгага Беагі ої тап а5 Бе Феїегтіпед ог 5осразйс ргосе55е5 15 єїмеп. 

Тре аїєогіфт ої дегегтіпацоп ої рагатеїег5я ої сагатас гру 15 Феуеіорей, Шаг 
ргоуідез Ше гесеїрі а5упаріойс еййсіепсу, азутріойс погплайіу апа 5ігопеЇу сопбізіепі 
езптапоп5 ої аггаує ої питап Беагі груп. Ваз5ед оп пефоа ої пахітаї ріайзіШіу, ап 
аїбогійпт 15 деуеіореа аг зиррозвійоп, Шаї сагадіо іпіегуа! е5итаїе сап Бе рге5епіеа, а8 ап 
ассідепіа! ргосе55 Ше ехресіеад уаїше ої у/рісі срапеез їп йте оп а рагтпопіоцз Іам, апда 
пої5е ра5 погта! ді5гібийоп. 

Ргеуміоця гезеагспез 5помеа Шаг аї ситепі е5ітайоп ої агауз ої сагаіас груп її 
15 епоцеб (0 ц5е Бе уагіабіез ої 5(апдага Іепеїб. Тіте ої апаЇузізє ої 5еРесйоп, уурісі Ба5 
іо 200 уаїце8, дое5 пої ехсеед 12 5ес. Ап аїєогійт апа ргоєгат ої деїегтіпацоп ої 
десгее ої Фегегтіпізт ої сагдіо група 15 Феусеіоред ас Фійегепі рзуспо рбузіоіовісаї 
егтлія раї аПом/5 вето Фе уаїше ої огвапігайоп ої Фе вої аггауз ої ритап реагі гпуфт. 
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